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C^ ■ Abstract 

Experimental data from the reaction of an 8.0 GeV/c 7r~ beam incident on a ^^^Au 
target have been analyzed in order to investigate the integrated breakup time scale 
for hot residues. Alpha-particle energy spectra and particle angular distributions 
supported by a momentum tensor analysis suggest that at large excitation energy, 
above 3-5 MeV/nucleon, light-charged particles are emitted prior to or at the same 
time as the emission of the heavy fragments. Comparison with the SMM and GEM- 
INI models is presented. A binary fission-like mechanism fits the experimental data 
at low excitation energies, but seems unable to reproduce the data at excitation 
energies above 3-5 MeV/nucleon. 
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During the past decade a major effort has been directed toward studies of 
the decay of hot nuclei. Experiments with energetic hght projectiles offer a 
unique perspective into decay processes at excitation energies up to about 
9 MeV/nucleon. The formation of highly-excited heavy residues in such re- 
actions is less complicated than in heavy-ion reactions, since only a single 
hot source is formed and small angular momenta and little compression are 
involved [1,2]. The subsequent evolution of the emission process for heavy frag- 
ments as a function of excitation energy is of major concern in understanding 
the decay mechanisms of hot nuclear matter. 

At low excitation energies fission fragments are emitted near the end of the 
sequential decay chain. At an excitation energy of about 3 MeV/nucleon the 
emission is delayed by about 10^ fm/c for the asymmetric split and this delay 
is much larger for symmetric division [3]. Recently it has been shown that 
even at E* /A of about 4 MeV/nucleon, fission fragments from the decay of 
hot gold-like nuclei are observed with the probability of about 20-40% [4]. The 
onset of different decay processes involving the prompt emission of multiple 
fragments is predicted to occur at excitation energies above 4-5 MeV/nucleon 
[5,6]. 

The experimental evidence for the transition from a sequential binary de- 
cay picture to a simultaneous breakup mechanism is discussed in this letter. 
Some other signals of the expected transition have recently been presented 
[7,8] on the basis of the experimental data set analyzed for the purpose of 
this letter. The Statistical Multifragmentation Model (SMM) [5,9] is one of 
the "explosive" models and the predictions from this model are compared with 
the experimental data. A freeze-out volume of V^ = 3- Vq (radius R = 1.44-_Ro)! 
without any extra expansion energy [8], which describes the Z > 6 multiplicity 
and charge distributions at large excitation energies, was assumed for all SMM 
calculations. Results from the sequential decay model GEMINI [10], includ- 
ing decay time scales and modified to calculate the trajectories in a mutual 
Coulomb field [11] are also presented. As input for all simulations, the mass 
and charge of the hot residual nucleus at a given excitation energy were based 
on event reconstruction procedure of the experimental data [12]. 

Experiment E900a was performed at the AGS Accelerator at Brookhaven 
National Laboratory. A secondary 8 GeV/c beam of tagged ir" irradiated 
a 2 mg/cm^ gold target. The ISiS Air detector array detects charged particles 
and fragments with charge up to Z f« 16 units and E/A si 1 — 92 MeV [13]. 
The experimental methods and previous results can be found in [12,14]. In 
this letter attention is concentrated on events that detect at least one heavy 
Intermediate Mass Fragment, IMF, specifically fragments with Zjmf = 8 — 16 
charge units. 

The probability for events with heavy IMFs is presented as a function of E*/A 



in Fig. 1. It is observed that the probabihty is very low at low E* /A, increases 
noticeably in the excitation energy range 3-6 MeV/nucleon and reaches about 
40% at about E* /A si 7 — 8 MeV/nucleon. Taking into account 74% detection 
efficiency, it can be concluded that the selected events are typical events at 
the excitation energy of about 6-8 MeV/nucleon, where the detected average 
number of thermal charged particles is about 12-15, the number of fast cascade 
particle is 10-15, based on a sample of 10^ events in this E* /A range. The 
thermal particles and fragments were selected by imposing a cutoff energy 
that accepted only fragments below 30 MeV for Z = 1 and 9 ■ Z + 40 MeV 
for heavier fragments [15]. SMM and GEMINI models predict for selected 
events that up to E*/A=7-9 MeV/nucleon the heavy fragment multiplicity 
with charge Z > 8 is two or slightly larger than two. The second, mainly 
heavier, fragment than that detected is too slow to be efficiently detected 
by ISiS. The event reconstruction procedure gives an estimation of the total 
missing charge that is consistent with a single missing heavy fragment which is 
equivalent to an IMF for E*/A > 7 MeV/nucleon. It is observed on Fig. 1 that 
experimental data are in good agreement with the SMM model predictions up 
to E*/A Ri 8 MeV/nucleon. The GEMINI model with a density parameter 
a=A/10 MeV~^ underestimates the probability to detect heavy IMFs above 
excitation energies 3-5 MeV/nucleon. 

In the 8.0 GeV/c tt^ -|-^^^ Au reaction the center-of-mass system, CM, is very 
close to the laboratory system, LAB, Vgource ~ 0.01c. It is observed that the 
thermal particles and fragments are emitted almost isotropically [16-18], con- 
sistent with emission from a thermal-like system. The data were also checked 
to insure that fast particles were not correlated with the heavy fragments and 
their spectra were independent of the angle relative to the heavy fragment 
axis. As shown in model calculations [1,2], the initial spin of the hot nuclei is 
low, so for the studies of the thermal decay of the hot nuclei the beam axis 
doesn't define any special direction. Therefore, the direction of the heaviest 
detected IMF in the LAB system was used as the reference axis for further 
studies. 

Alpha-particle kinetic energy spectra are presented in Fig. 2, for a few se- 
lected angles, Qrei, relative to the axis of the heaviest detected fragment 
ZiMF = 8—16. On the left side of Fig. 2 the spectra at E*/A=l-3 MeV/nucleon 
are presented. Over most of the angular range the shape of the alpha-particle 
energy spectra doesn't change much. However, the thermal alpha-particle aver- 
age kinetic energy, < E^in >, has a maximum value at small- and large- angles, 
< Ekin >~ 23 MeV. The minimum is observed at 50 deg, < E^in >~ 20 MeV. 
The GEMINI model and SMM model (running at E*/A < 3 MeV/nucleon 
essentially as a sequential evaporation model) ffi the shape of the spectra, 
however, predict a more pronounced evolution than the experimental data. 

The spectra at E*/A=5-7 MeV/nucleon are presented on the right side of 



the Fig. 2. It is important to note that the multifragmentation model SMM 
reproduces the slope of the spectra up to about 60 MeV. The GEMINI model 
predicts much steeper spectra. It is also observed that at angles close to the 
axis defined by a heavy fragment, the alpha-particle kinetic energy spectra 
have a different shape than the spectra at the larger angles. The average ki- 
netic energy is < Ekin >~ 29 MeV, compared to a nearly constant value of 
< Ekin >~ 24 MeV at the angles between 50 and 170 deg. Similar evolu- 
tion of the energy spectra is also observed for protons and lithium isotopes, 
i.e. at large excitation and small angles, Orel, the kinetic energy spectra are 
shifted towards higher energy. The modified GEMINI model [11] predicts a 
shift similar to the observed one. Tracing the calculations, it was found that at 
E* /A=5-7 MeV/nucleon GEMINI predicts a short, exponentially distributed 
time for heavy fragment separation, on average about 100 fm/c. According to 
the model, most of the alpha-particles are emitted from the hot fragments. The 
modified version of GEMINI model [11] shows that preferentially low energy 
alpha-particles are suppressed towards small angles. Ore/- Particles emitted 
from the complementary fragment towards the selected fragment are defiected 
by the Coulomb field to larger angles Qrei- This result of the model indicates 
that secondary decay of the hot fragments (emitted promptly after thermal- 
ization) is able to explain the kinetic energy shift at small angles, in despite 
of all doubts concerning the validity of GEMINI model at large excitation 
energies. 

The SMM model doesn't show the observed energy shift of alpha-particles at 
small angles. This failure seems to be related to the difficulties to define a set 
of parameters at the freeze-out point configuration. In this context it should 
be noticed that the version of the SMM model used assumes that secondary 
decay occurs after the heavy fragments have been fully accelerated, at infinite 
separation distance, taking into account only two-body Coulomb interaction. 

It is interesting to discuss the possibility to describe the shift of alpha-particle 
kinetic energy at large excitation energy, and at small angles, Orel, by a fission- 
like decay scenario. This requires that the hot nucleus survives as a single 
self-bound object, the heavy fragments separate after a certain delay time, 
and the alpha-particles are mainly emitted prior to fragment separation, pos- 
sibly from a strongly deformed system. One can envision two touching hot 
spherical nuclei at rest, which would resemble a largely deformed system. In 
a such a case the alpha-particle energy depends only slightly on the emission 
angle. Particle emission from the neck at the scission configuration seems to 
be the only mechanism to enhance the yield of low energy particles emitted 
perpendicular to the fragments. However, it is observed that the average emit- 
ted thermal alpha-particle multiplicity at E*/A=5-7 MeV/nucleon is about 5 
and to suppress low energy alpha-particles at small angles, the emission from 
the neck region would have to be a dominant decay process. In other words, 
the hot nucleus should emit multiple particles mainly from the strongly de- 



formed system, primarily from the neck region, and not from the tip region. 
Hence, it seems to be difficult to describe such behavior by a fission-like decay 
mechanism. 

The angular multiplicity distributions of thermal particles and light fragments 
{Z = 1 — 7) relative to the emission axis of a Zjmf = 8 — 16 fragment are shown 
on the left side of Fig. 3. The particle multiplicity distributions are plotted for 
selected excitation energy bins, along with SMM model and GEMINI model 
predictions. In all such angular distribution plots the decrease at angles below 
15 deg refiects the ISiS angular granularity. The GEMINI model and SMM 
model at E*/A below 3 MeV/nucleon predict a bump at small angles due to 
particles emitted from the accelerated fragments. The absence of enhanced 
emission along the heavy IMF axis in the experimental data indicates that 
the accelerated fragments have lower excitation energies than predicted by 
the models. Such a picture is consistent with the binary fission scenario, i.e. 
the fragments are emitted at the end of the sequential decay chain, in a good 
agreement with refs. [3,4]. 

The angular distributions at larger E* /A show that the probability of detect- 
ing a particle or fragment close to the selected fragment is suppressed, and 
suggests that most of the charged particles are emitted almost simultaneous 
with the selected fragment emission. In other words, the Coulomb field of 
the heavy fragment makes a cone-shadow not accessible to the other particles 
emitted at the same time. Consistent with this picture, the SMM model predic- 
tions filtered with the detection efficiency procedure describe the experimental 
data well. At the very large excitation energy, E* /A=7-9 MeV/nucleon, the 
GEMINI model also fits the data due to the very short time scale predicted by 
the model: the heavy fragment relative separation time of about 75 fm/c, and 
the first neutron emission decay time of about 6 fm/c. However, the GEMINI 
predictions should be taken with the special care at large excitation energies. 
The predicted decay time at E* /A=l-9 MeV/nucleon is shorter than the time 
to pass the distance equal to the nucleus diameter with a speed of light. This 
indicates that the decay time predicted by the GEMINI model is shorter than 
the energy relaxation time, which contradicts to the model assumption. 

In order to investigate systematically the total breakup time of the system 
and the hypothesis of simultaneous decay, a momentum tensor analysis [19] 
was performed on the charged particles. The tensor was built for each event 
from all detected thermal charged particles and light fragments with Z = 1 — 7 
charge units. Each element of the tensor, Qij, was calculated as: 



n — \^ i^j__t^j__ 



where, N is the total number of particles in an event, p^"^ is the particle mo- 
mentum in the center-of-mass system of A^ particles with Z = 1 — 7 and p\j 
are the particle momentum Cartesian coordinates. The tensor was diagonal- 
ized and the angle between the vector associated with the largest eigenvalue 
and the emission direction of the selected heaviest fragment, Zjmf = 8 — 16, 
was defined as momentum flow angle, Qfiow This method permits study of the 
correlations between the heavy fragments and lighter products of the reaction. 

The distributions of the angle Qfiow are presented on the right side of the 
Fig. 3. The experimental data are almost isotropically distributed at low ex- 
citation energy but develop a maximum at 90 deg that becomes increasingly 
pronounced with increasing excitation energy. An isotropic Qfiow distribution 
assumption describes more than 95% of events at E* /A=l-3 MeV/nucleon 
but less than 60% of events at E*/A=7-9 MeV/nucleon. These experimental 
results should be discussed in a frame of different decay scenarios. 

If particles are emitted isotropically, before fragment separation, then Qfiow 
has also an isotropic distribution. According to the binary flssion scenario, this 
should be the case for the experimental data at E*/A below 3 MeV/nucleon. 
The small enhancement of events, below 5%, with Qjiow greater than zero 
seems to be due to particle emission perpendicular to the flssion axis, perhaps 
from the neck region at the scission conflguration. 

If there are additional particles emitted from the accelerated fragments, then 
the momentum ellipsoid should be elongated in the direction of the heavy 
fragment emission and the Qfiow distribution should be focused at small angles 
(the kinematical focusing effect). This behavior is shown for low excitation 
energy by the GEMINI and SMM predictions, corresponding to the angular 
distributions on the left side of the Fig. 3. 

On the other hand, if the associated charged particles are emitted at the same 
time as the heavy fragment, then preferential Coulomb focusing in a direction 
perpendicular to the selected heavy fragment emission axis is expected. In such 
a case the momentum ellipsoid is elongated perpendicular to the heavy frag- 
ment emission direction and Q/iow should present a maximum at 90 deg. This 
focusing effect normal to the heavy fragment axis gives evidence that at high 
excitation energy the Coulomb focusing effect is stronger than any kinematical 
focusing effects. The prompt decay scenario, including the Coulomb focusing 
effect, is contained in the SMM model. Thus, at large excitation energy the 
SMM model describes the Qfiow distribution well. The modifled version of the 
GEMINI model, with very short emission times after thermalization also flts 
the the Qfiow data at large excitation energy. 

The integrated time scale for breakup of the thermal-like system at excitation 
energies above 5-7 MeV/nucleon can be estimated from the experimental data. 



At high excitation energy the observed stronger Coulomb focusing relative 
to kinematical focusing (right side of the Fig. 3) means that fewer particles 
are emitted after the acceleration time compared to the number of particles 
emitted within the first 200 fm/c after fragment formation. The modified 
GEMINI model gives an estimation that more than about 75% of particles 
emitted after fragment separation would have to be emitted within the first 
200 fm/c to create stronger Coulomb focusing effect than kinematical focusing 
effect. The light fragment, Z < 7 , angular distributions (Fig. 3) and alpha- 
particle kinetic energy spectra (Fig. 2) indicate that not many particles can 
be emitted after the thermalization but before heavy fragment separation. 
Taking into account these experimental results, it can be concluded that at 
excitation energies above 5-7 MeV/nucleon the disintegration of the system 
typically takes place within a total time of approximatively 200 fm/c after 
thermalization. This short lifetime is in a good agreement with the results 
obtained from the IMF-IMF intensity-interferometry analysis based on the 
same experimental data [7]. 

The presented results are obtained for events with detected heavy IMF, Z = 
8 — 16, and IMF-IMF interferometry analysis was based on the set of events 
with detected two light IMFs, Zi^2 = 4 — 9. At a large excitation energies, 
above 5-7 MeV/nucleon these two set of events represent a large, almost com- 
plementary set of events, and cover together more than 70% of detected events. 
However, it seems that the method based on the sample of events with the 
splits closer to a symmetric split is a better choice looking for the experimental 
evidences of the expected transformation from the shape deformation decay 
mode to the prompt decay mode. Moreover, the set of events with detected 
heavy fragments covers a large range of excitation energy, the fission-like events 
were already successfully used in the past [3,20] to study the decay time scale 
at low excitation energy, and the IMF-IMF interferometry technique has ac- 
cess only to a marginally low fraction of events at low excitation energy due 
to a very low IMF emission probability. 

In conclusion, it is shown that for hot nuclei at excitation energies below 3-5 
MeV/nucleon the experimental data are consistent with a sequential binary 
breakup scenario [3,4], i.e. the heavy fragments are emitted at the end of a 
sequential decay chain. The fission-like decay scenario above E* /A of about 
6 MeV/nucleon seems unable to fit the experimental data, because it would 
require emission of too many particles from the neck of the scission configura- 
tion. At excitation energies above 3-5 MeV/nucleon, charged-particle emission 
appears to occur nearly simultaneously with the emission of heavy fragments, 
ZiMF = 8 — 16. The integrated breakup time of the system estimated for the 
experimental data is approximatively 200 fm/c. This observation is in agree- 
ment with the delay times from IMF-IMF small angle correlation studies [7] 
and Statistical Multifragmentation Model (SMM) predictions. 



Acknowledgements: This work was supported by the U.S. Department of 
Energy and National Science Foundation, the National and Engineering Re- 
search Council of Canada, Grant No. P03B 048 15 of the Polish State Com- 
mittee for Scientific Research, Indiana University Office of Research and the 
University Graduate School, Simon Eraser University and the Robert A. Welch 
Eoundation. 



References 

[I] J. Cugnon, Nucl. Phys. A 462 (1987) 751, and J. Cugnon et al, Nucl. Phys. 
A 470 (1987) 558. 

[2] A.S.Botvina et al, Phys. Lett. B 205 (1988) 421. 

[3] D. Hilscher et al, Phys. Rev. Lett. 62 (1989) 1099. 

[4] U. Jahnke et al, Phys. Rev. Lett. 83 (1999) 4959, and L. Pienkowski et al, 
Nucl. Phys. A 655 (1999) 269c. 

[5] J. P. Bondorf, A.S. Botvina, A.S. Iljinov, LN. Mishustin and K. Sneppen, Phys. 
Rep. 257 (1995) 133. 

[6] D.H.E. Gross, Rep. Prog. Phys. 53 (1990) 605, and D.H.E. Gross, Phys. Rep. 
279 (1997) 119. 

[7] L. Beaulieu et al, Phys. Rev. Lett., 84 (2000) 5971. 

[8] T. Lefort et al, Phys. Rev. C 62 (2000) 31604. 

[9] A.S. Botvina, A.S. Iljinov and I.N. Mishustin, Nucl. Phys. A 507 (1990) 649. 

[10] R.J. Charity et al, Nucl. Phys. A 483 (1988) 371, GEMINLcode obtained via 
ftp from wunmr.wustl.edu. 

[II] W. Gawlikowicz, Acta Phys. Pol. B 28 (1997) 1687. 
[12] L. Beaulieu et al, Phys. Lett. B 463 (1999) 159. 

[13] K. Kwiatkowski et al, Nucl. Instr. Meth. A 360 (1995) 571. 

[14] T. Lefort et al, Phys. Rev. Lett. 83 (1999) 4033. 

[15] K. Kwiatkowski et al, Phys. Lett. B 423 (1998) 21. 

[16] T. Lefort et al. Proceedings of the XXXVII International Winter Meeting on 
Nuclear Physics, edited by I.Ori Bormio (Italy) 25-30 January 1999, page 7. 

[17] F. Goldenbaum et al, Phys. Rev. Lett. 77 (1996) 1230. 

[18] W.-c. Hsi et al, Phys. Rev. C 60 (1999) 034609. 

[19] J. Cugnon and D. L'Hote, Nucl. Phys. A 397 (1983) 519. 

[20] D. Hilscher and H. Rossner, Ann. Phys. Fr. 17 (1992) 471. 




8 10 

E*/A (MeV) 



Fig. 1. Dots present the probability to detect an event with a heavy IMF, 
ZiMF = 8 — 16 and sohd (dashed) hne presents the SMM (GEMINI) model predic- 
tions filtered with the experimental detection efficiency. An initial angular momen- 
tum of hot nucleus L = 20/i was assumed for GEMINI model calculations. 
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Fig. 2. Alpha-particle energy spectra for selected excitation energy and angular 
bins relative to the heaviest fragment, Zjmf = 8 — 16, emission axis, Qrei- Solid 
(dashed) lines present the SMM (GEMINI) prediction filtered with the experimental 
detection efficiency. The simulation results are normalized to the experimental data 
assuming the same cross section at a given angle. 
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Fig. 3. Laboratory angular distributions with respect to the axis defined by the emis- 
sion of the heaviest detected IMF, Zjmf = 8 — 16 charge units, for the selected ex- 
citation energy bins. Left side: Average multiplicity of the detected thermal charged 
particles with Z = 1 — 7 charge units as a function of OreZ- Right side: probability 
distribution of the momentum flow angle, Qfiow- Solid dots present the experimental 
data. Solid (dotted) lines present the SMM model predictions filtered (non filtered) 
with the experimental detection efficiency. Dashed lines are from GEMINI model 
predictions filtered with the experimental detection efficiency. 
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